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ABSTRACT 


Since regular operation of the DOE/NASA Mod-1 wind turbine began in October 1979 about 10 nearby 
households have complained of noise from the machine. Ctevelopnent of the NASA-LeRC wind turbine sound 
prediction code began in May 1980 as part of an effort to understand and reduce the noise generated by 
ffed-I, Tone sound levels predicted with this code are in generally good agreement with measured data 
taken in the vicinity f^d-1 wind turbine (less than 2 rotor diameters). Comparison in the far field 
indicates that propagation effects due to terrain and atmospheric conditions may be amplifying the actual 
sound levels by about 6 cto. Parametric analysis using the code has shown that the predominant contri- 
butor to Mod-1 rotor noise are (1) the velocity deficit in the wake of the support tower, (2) the high 
rotor speed, and (3) off-optimum operation. 


INTRCXXJCTION 


'Since' regular operation of the Mod- 1 wind 
turbine began in October 1979 about ten 
households have complained of noise from the 
machine. The character of the noise is 
descri bed as^an audible ”thum p” at a repetition 
rate equal to the blade passing frequency. In 
some instances, low frequency acoustic energy 
has resulted in complaints of vibrations within 
the homes. 


Since January 1980, efforts have been directed 
at identifying the causes of the noise ar«j 
methods of reducing it. One effort has been the 
development of computer prediction codes and 
obtaining experimental data to verify these 
codes. The objective of this effort is to 
understand the noise generated by Ntod-1 and to 
prevent it from being a problem on advanced 
machines. Development of a wind turbine sound 
prediction code began at the NASA-LeRC in May 
198G. This report presents (1) measured data 
that characterize the noise problem at Mod-1, 

(2) the analytical method used in the NASA LeRC 
code, (3) validation of the code using experi- 
mental data from Ntod-1, and (4) applications of 
the code tc the Mod-1 situation, 

SUMMARY 

Measured data taken at the Mod-1 site show that 
the impulsive character of tne noise is composed 
of harmonics of the blade passing frequency. 

While some of these harmonics exist below the 
audible frequency range (about 20 to 20,000 Hz), 
many are aoove the nominal audible threshold of 
20 Hz. 

A computer program (WTSOIM)) has been developed 
for calculating the intensity and frequency 
characteristics of sounds generated by wind 
turbines in a non-uniform wind flow field. The 
results calculated with this code are generally 
in good agreement with Mod-1 measured sound spec- 
tra. However, propagation effects due to terraiir 
-and atmos^eric' cronditions have complicated the 
amplitude correlation with l^d-1 data in the far 
field. These effects have been estimated by the 
code to cause an amplification of 6 dB or more 
at a home that has registered complaints. 


The WTSOUND code has been used to determine the 
source of the noise generation from Mod-1 as 
well as Identifying operating conditions 
associated with the highest noise levels, 
hfejdeling of the wind flow field characteristics 
shows that the predominant source of noise from 
Nfod-1 is the wind velocity deficit in the wake 
of the tower. Because the rotor plane is 
downwind of the tower, this deficit produces 
charges in the aerodynamic forces on the blades 
resulting in sound pressure variations in the 
acoustic field. The level of the sound pressure 
variations is mst strongly a function of rotor 
speed and wind speed. Reducing the rotor speed 
from 35 to 23 rpm is predicted to reduce sound 
levels by about 11 d8. The increase in sound 
levels with windspeed is expected to be about 12 
dB between cut in and rated windspeeds. Vari- 
ation in sound level with power is not clear-cut 
since hi^ sound levels can occur during both 
high and low power conditions. 

MOO-1 MEASURED DATA 

In January 1980, the Solar Energy Research 
Institute (SERI) and the General Electric 
Company under DOE funding began visiting the 
Ntod-1 site to investigate the noise problem. 
During those visits, tape recordings of the 
sound pressures were made both near the Mod-1 
and at residences that had complained of noise. 

A good example of the recorded time history of 
the pressure variation near the Mod-1 recorded 
by SERI is shown in Figure 1. The impulsive 
nature is clearly apparent as each blade passes 
behind the tower structure. The passage of the 
blade throu^ the wake of each leg of the tower 
is seen as the two upper peaks approximately 
0.075 seconds apart. 

To. bet ter understand the character of the 
pressure impulse, spectrum analyses were 
performed on the measured data. One sample 
.spectrum of the sound near the Mod-1 is shown in 
Figure 2. This was recorded when the machine 
was generating a heavy ”thumping" sound at a 
residence. From this analysis, it is obvious 
that the impulsive character is composed of 
harmonics of the blade passage frequency (1,16 
Hz), Che hipest amplitude harmonics 
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are in the sub-audible frequency range (less 
than about 20 Hz), tones are apparent in the 
audible frequency range above 20 Hz. 



Figure 1 - Time history of the sound 
pressure variation near 
the Mod-1 



FREQUENCY, HZ 

Figure 2 - Spectrum Analysis of 
sound variation near 
the Mod-1 


ANALYTICAL METHOD OF THE SOUND CODE 


wind turbines in non-uniform flow fields. 

The predominant sound produced by a wind turbine 
is associated directly with the aerodynamic 
pressures on the blades. These pressures can be 
related for convenience to the thrust and torque 
forces on the rotor. The thrust and torque 
forces have components that are both steady and 
unsteady in time. The steady forces produce 
sound called rotational noise, which consists of 
pressure variations in the acoustic field at the 
blade passing frequency with harmonics of 
rapidly decreasing magnitude. The unsteady 
forces may be either periodic (i.e. tower shadow 
and windshear) or random (i.e. gusts). Noise 
due to periodic unsteady forces may be dominant 
over rotational noise and generated higher 
harmonics of amplitude comparable to that of the 
fundamental. 

The method used to determine the sound pressure 
levels in the acoustic field is described by the 
flow chart in Figure 3. This procedure can be 
summarized as follows: (1) calculation of the 
steady aerodynamic blade forces, (2) variation 
in these forces due to unsteady aerodynamics, 

(3) Fourier analysis of the force variation, and 
(A) calculation of sound pressure levels in the 
acoustic field. These steps are explained in 
more detail in the following sections. 


INPUT CALCULATION OUTPUT 



The approach* used in developing the wind turbine 
sound prediction code, WTSOIND, was to apply an 
avail^le theory used for calculating noise from 
conventional aircraft propellers. The 
development of such theories goes back to 1937, 
when Gutin first successfully calculated the 
noise from a propeller in a uniform flow field 
(ref. 1). Since then, Gut in ^s theory has been 
extended to include the effects of non-uniform 
flow fields (ref. 2) and applied to helicopters 
and turbomachinery as well as propellers. The 
WT5CUND code was written using this theory to 
provide a means of calculating sound intensity 
and frequency characteristics specifically for 


Figure 3 - Flow chart of sound 

pressure level calculation 


Steady Aerodynamic Blade Forces 

The total thrust force and torque on a rotor in 
uniform flow is determined from blade element- 
momentum theory. The reader is referred to 
reference 3 for development of this theory. 
There are several computer programs available 
that use blade element- momentum theory, 
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(3a) 


including the PROP code (ref. 3). The PROP code 
requires modeling of the wind turbine 
characteristics and operating conditions (i.e. 
planform, twist, rpm, windspeed, etc) to 
calculate the steady torque and thrust force on 
the rotor. The total torque and thrust on the 
rotor can then be resolved into equivalent 
forces actir^ on each blade at a single point. 
The point used in this analysis is the 75 
percent’ radius of the blade. That radial 
location is also used when modeling the unsteady 
■|erodyr^mic variation. 

Unsteady Aerodynamic Blade Forces 


Crce the steady forces have been determined, the 
unsteady forces are obtained throu^ perturb- 
ation with non-dimensional force coefficients. 
Figure 4 shows the the wind velocity and force 
vector dragram of an airfoxl element of chord 
length c . This airfoil is operating at a pitch 
angl^^B with respect to the plane of rotation. 
The velocity of the wind at the rotor plane, , 
confines with velocity due to rotation, Vq, 
to give the relative velocity vector, Vpt. The 
relative velocity vector acts at an angle of 
attack ct with respect to the airfoil chord 
line. The relative velocity and angle of attack 
are given by 



The airfoil lift and drag coefficients are shown 
as the heavy line vectors in Figure 4. 



Figure 4 - Wind velocity and force 
vector diagram 


These coefficients can be transformed into a 
thrust force coefficient, Cj , acting perpendi- 
cular to the rotor plane and a torque force 
coefficient, Cq , acting parallel to the rotor 
plane where 


cos (a + 9) -f Cjj sin(a + 9) 


and 

Cq * sin(a + 9) - cos(a + 6) 


(3b) 


me thrust and torque coefficients can now be 
used to determine the unsteady forces associated 
with periodic variations in the wind velocity, 

. Such a periodic variation will occur as 
the blade rotates through wind shear or behind 
the wind turbine support tower. The quasi- 
steady state values of blade thrust force, T^ , 
and torque force, , at any rotor azimuth 
position, tJj , are given by 


and 




c'*' 

Q s 


(4a) 


(4b) 


in which the superscript s denotes the steady 
forces and coefficients respectively. 

Fourier Analysis of the Force Variation 

The next step is to a perfqrn^ a Fourier analysis 
of the blade force variation. The complex 
Fourier coefficients for the thrust and torque 
forces respectively are given by 

/ 2t/Q ^ , 

eiP^V(T)dT (5a) 


/•2-n/Q 

•J-n / (56) 

in which n is the rotor speed. 

These coeffi cients are determined in the program 
using the IBM siiD routine FORIT. FORIT gives 
real Fourier coefficients Ap and Bp for the 
cosine and sine terms respectively. The 
transformation to complex coefficients is 

a * A + iB^ (6a) 

P P P 

or 

a =» A - iB , ^ 

-P P P (6b) 

for 

p _> 0 

A correction to the quasi-steady state analysis 
can now be made by including the effects of 
LTisteady aerodynamics. The approach used here 
to determine the response of the airfoil was 
developed by Sears (ref. 4). The correction is 
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VALIDATION OF THE WT50UND CODE 


given by a simple expression called a Sears 
function, which is used as a factor to the 
Fourier coefficients. The expression used to 
approximate the Sears function is 

S<a) . - 2(1 ^2.,)]} 

(1 + 2tt0)1/2 ( 


for 


0 > 


0 


where 


0 


pOc 

2V 


(7b) 


The coefficients an and a_p in Equations ( 6 ) 
are multiplied by the Sears factor before use in 
calculating sound pressure levels. 

Sound Pressure Levels in the Acoustic Field 


The mathematical relationship for calculating 
the sound pressure levels from the Fourier 
coefficients of the blade force variation was 
obtained by Lowson (ref, 5), The RMS pressure 
variation of the harmonic of the blade 
passage frequency is given by the following 
equations: 


P 

n 


Atts 



ip(4^“TT/2) 

nB-p 


Ck 


sin y) 


To verify the accuracy of the WTSOUND code, 
measured data from the Mod-1 wind turbine were 
used. The data presented here were taken from 
references 6 and 7, 


Correlation with the Mod-1 Measured Data 

On June 10, 1980 at 12:36 a.m, sound levels were 
being measured with a microphone located about 
240 ft from the Mod-1. The wind turbine was 
operating at 34.6 rpm in a 30 mph wind and 
generating about 750 kw Into a load bank. The 
measured sound spectrum is shown in Figure 5 and 
has been designated case GE 180. These measured 
sourxj levels consist of both tones and broad- 
band noise. The tone levels, shown as the 
narrow peaks, “are generated by the wind 
turbine. The broad-band level, shown as the 
flat areas of the spectrum, is composed of 
ambient wind noise as well as blade vortex noise. 





cos y - 


(8a) 


lp(((.-7r/2) 


X |a'^ cos Y - P 

“P -P 


^ (t® cos Y Q® 


and ( 8 b) 

n to 

in which 

B is the number of blades 

S is the distance from the rotor 

Y,(|) are azimuth and altitude angles 
to the listener 

r is the blade radius where the thrust 
and torque forces are assumed to act 


FREQUENCY, HZ 


Figure 5 - Mod-1 measured sound 
levels for case 6 E 180 


Since the WTSOUND code calculates tone sound 
pressure levels, the effect of the broad-band 
noise on the spectrum must be removed for 
comparison with the code. To do this, two 
points on the spectrum were chosen that were 
believed to be at the broad-band noise level 
( 8 i and 62 in Figure 5), The broad-band 
levels throughout the spectrum were calculated 
using the following empirical model: 



where 

L^ is the broadband sound level in 
' cB at frequency f 

L^Lp 3 re the broadband sound levels 
^ ^ in dB at points Bx and B 2 


J is the standard Bessel function 
Co is the speed of sound 


tl, f 2 are the frequencies corresponding 
to points Bx and B 2 
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The tone-only levels were then determined by 


= 10 log 110 ^ - 10 ^ I (iO) 


where 


is the tone sound level in cB at 
: : frequency f 

: :: is the measured total sound level in 

^ dB at frequency f 

Some typical adjusted tone-only sound pressure 
levels are shown as the squares in Figure 5. 
fstote that this correction is small and can only 
be seen in areas where the tone and broad-band 
levels are comparable. 

An analytical model of the Mod-1 wind turbine 
was de veloped for case G£ 180. Steady thrust 
"and torque forces, required input for the 
WTSOUND code, were calculated using the PROP 
p^rTbf marc e code. The wind velocity deficit in 
of the tower was approximated as an 
average cf the velocity profiles at the 69, 75, 
and 81 percent blade radius. The velocity 
profile w^e t^en from scale model wind tunnel 
tests of the Mod-1 tower (ref, 8). Figure 6 
shows the neasured wake velocities as well as 
the analytical wake model assumed. Note that 
the measured and assumed profiles were taken at 
different distances downwind of the tower. No 
correction to the wake for distance appears to 
be rej^uired, and none was made. 


lines are the tone levels predicted by the 
WTSCXih® code. In general the_code predicts the 
amplitudes of the hipest harmonics very well. 
Also the roll-off rates (arrplitude decreasing 
with frequency) of the harmonics compare 
favorably . The differences _.that do exist are 
believed to be associated with the tower shadow 
model mentioned earlier. 



FREQUENCY, HZ 

Figure 7 - Comparison of theoretical 
and measured sound spectra 
for case GE 180 



BLADE AZIMUTHAL ANGLE, DEG. 
Figure 6 - Mod-1 tower shadow 


Lowscn^s equation gives free-space sound 
pressures with no effect of reflection from 
nearby solid bodies. A 6 dB increase has been 
included in the WT5CUND code analysis to account 
for reflection when using microphones near the 
grourd. This correction is common practice when 
calculating airplane propeller noise at ground 
level (ref. 9). 

Figure 7 shows the comparison of neasured and 
predicted sound levels. The. squares are the 
same as those shown in Figure 5. The vertical 


To more conveniently characterize the overall 
sound level of this spectrum, the root sum 
square (RSS) of the harmonics between 20 and 50 
Hz was calculated. The lower end of this range 
was chosen as the nominal threshold of hearing. 
The upper bound was chosen because the tone 
levels above this frequency fall below the 
broad-bana noise. As shown in Figure 7 the 
20-50 Hz RSS sound pressure level given by the 
WTSOUND code is 2 dB lower than the measured 
data. This agreement is felt to be very good. 

Table 1 summarizes the results from GE 180 as 
well as two other examples, ca^es GE 140 and GE 
IB. These two cases were chosen to test the 
ability of the code to predict changes in rotor 
speed and distance from the machine, 

GE 140 documents the sound levels near the Mod-1 
on June 9, 1980 at 9:22 p.m. The wind turbine 
was operating at a reduced speed of 22.7 rpm 
while again generating 750 kw. The wind speed 
of 26 mph was lower than the 30 mph wind of GE 
180. The measured RSS sound level for QE 140 is 
17 cB less than GE 180. Note that this 
reduction is exactly the same as predicted by 
the code. 

G£ IB documents the sound levels at some 
distance from the Mod-1 (3100 ft. downwind) on 
March 31 at 12:26 pm. This location is near a 
residence that has complained of noise from the 


machine. At this time the Mod-1 was operating 
at 34. 7 rpm in a 30 mph wind and producing about 
1850 kw into the utility grid. The measured 
20-50 Hz RS5 sound level is 21 dB lower than GE 
180. The predicted level, however, was 27 dB 
lower. Thus the actual level is about 6 dB 
hi^er than predicted by normal spherical 
dispersion. This 6 dB increase is believed to 
be associated with focusing of the sound due to 
terrain and atmospheric conditions. Though the 
amount of amplification may vary, this case 
demonstrates the use of the code to quantify the 
effects on sound propagation. 


Table I - Summary of Mod-1 
sound data cases 


ffiAR FIELD (1,3 D) 


CP&E 

m 

WIND SPEED. 

20-50 HZ RHS SOIWD PRESSURE LEVEL. DB | 

RPH 

CODE 

MEASURED 

DIFFERENCE 

GE 180 

Vi.b 

30 

79 

81 

-2 

GE IMO 

22 J 

26 

62 

69 

-2 

m FIELD (15.5 D) 

GE IB 

M.7 

30 

52 

60 

-8 


APPLICATION OF ThE SOUND COCX 

After gaining confidence in the ability of 
WTSOUND to predict the sound characteristics of 
wind turbines, the code was used to investigate 
the Mod-1 problem, in particular, studies were 
made of (1) the contribution to the sound level 
by characteristics of the flow field, and (2) 
the operating conditions that affect the Mod-1 
sound problem. 

Effect of Flow Field Characteristics 

The WTSOUND code was used to determine the 
relative contribution to the sound level by the 
characteristics of the flow through the rotor. 
Figure 8 shows the effect of the uniform flow, 
wind shear, and tower shadow. As stated 
earlier, the uniform flow field contribute 
significantly only to the fundamental blade 
passing frequency. Similarly, wind shear, which 
is essentially a 1/rev. variation, contributes 
only to the first few harrnonics. However, with 
the addition of tower shadow, many harmonics are 
introduced that have hi^ amplitudes up into the 
audible frequency range. Thus, the predominant 
contribution to the noise problem of Mod-1 is 
the wind velocity deficit in the wake of the 
tower. 



FREQUENCY, HZ 

Figure 8 - Effect of wind shear and 

tower shadow on wind turbine 
sound level (Mod-1 case GE IB) 

Effect of Operating Conditions 

In an effort to better understand the conditions 
under which the highest sound levels are 
produced, the WTSOUND code was used to study the 
effects of various operating parameters. In 
particular, the effects of power, rotor speed, 
wind speed, and direction were investigated. 

Figure 9 shows the predicted sound levels versus 
power output for both 35 and 23 rpm operation. 

The solid line shows normal operation as the 
wind speed varies from cut in through rated to 
cutout wind speed (6, 15, and 20 m/sec 
respectively). This analysis shows about a 12 
dB increase in sound level from cutin to rated 
wind speeds for both 35 and 23 rpm operation. 

This agrees reasonably well with empirical 
modeling of the measured data (ref. 7) that 
indicates a variation of about 14 dB over those 
wind speeds. 

During normal operation at 23 rpm, the sound 
levels predicted in Figure 9 are about 11 dB 
lower than those for 35 rpm. This also agrees 
well with measured data taken in June 1980. 

During those tests, the Mod-1 was operated at 
various speeds into a load bank. Analysis of 
these data showed a reduction of about 10 dB at 
a rotor speed of 23 rpm (ref. 7). 

The shaded area in Figure 9 represents operation 
below the available power in the wind for wind 
speeds less than rated. Operation in this 
region is characterized as off-optimum. Under 
this situation the pitch of the blade is 
increased toward feather, to limit power below 
the rating of the machine. The shaded area 
shows that hi^ sound levels may be produced 
even at low power levels. This conclusion seems 
to be supported by the fact that GE 180 (750 kw, 

30 mph) is further from normal operation than GE 
140 (750 kw, 26 mph) and is higher in sound 


416 






i 

1 



level than the 11 cB difference associated with 
rpm only« 

The sound levels at off-optimum operation above 
r ated have not been determined. Irrprovements in 
the_ unsteady aerodynamic force calculation in 
the WTSCXJNO code are being incorporated to 
analyze this area of operation. 



POWER, m 

Figure 9 - Effect of power, rpm, and 
windspeed on Mod-1 wind 
turbine sound levels 


The WTSOUHD code can also be used to calculate 
the directivity pattern of the sound. Figure 10 
shows the directivity pattern in a plane 1100 
ft. belo w the Mod-I hub height. The direction 
of highest sound pressure is directly downwind; 
however7^t~Is~ nearly as high upwind, Minirwim 
sound generation occurs in the rotor plane where 
levels are about 18 cB lower than downwind. 

This compares reasonably well with measured data 
in the near field showing about 15 cB variation 
(ref, 7). 


CONCLUSIONS 

The WTSOUND computer code shows generally good 
agreement with sound spectra measured in the 
vicinity of a wind turbine. In the far field, 
however, correlation of the absolute amplitude 
of the sound level is complicated by propagation 
effects. For the case in this study, terrain 
and meterological conditions caused an increase 
of about 6 dB. 

Analysis using the SOUND code shows that the 
predominant contributor to the noise problem of 
is the wind velocity deficit in the wake 
of the tower. Changes in the aerodynamic 
forces, as the blades pass throu^^ the deficit, 
produce sound pressure variations in the 
acoustic field. 

The level of the sound pressure variations are 
most directly affected by rotor speed and 
windspeed. Reducing the rotor speed from 35 to 
23 rpm is predicted to reduce sound levels by 


about 11 cB. The increase in sound levels with 
windspeed is predicted to be 12 dB between cutin 
and rated. 

ROTOR PLANE 

K 


I 

f 



Figure 10 - predicted sound directivity 
pattern In a plane 1100 ft. 
below the Mod-1 
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QUESTIONS AND ANSWERS 
L.A. Viterna 


From: G. Greene 

Q: What was the averaging time of the data you used for comparison? 

A: Approximately 5 minutes , 

From: D,W. Thomson 

Q: Are you working on a code for f < 20 Hz? The observational data seemed to strongly 

support structural excitation as a major source of annoyance. 

A: The code predicts harmonics starting at the blade passage frequency 1 Hz). We 

are using the 20-50 Hz RSS sound level to characterize audible annoyance . 

From; N.D. Kelley 

Q: Why does your model predict so many spectral peaks between blade passage and 50 Hz 

compared to OBS? 

A: This is due to the actual tower wake shape being different than the assumed wake 

based on scale model wind tunnel tests. We feel however^ that the assumed wake 
is adequate to acceptably predict the spectrum characteristics . 
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